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Abstract 



A very natural construction of integrable extensions of soliton systems is presented. The 
extension is made on the level of evolution equations by a modification of the algebra of 
dynamical fields. The paper is motivated by recent works of Wen-Xiu Ma et al. (Comp. 
Math. Appl. 60 (2010) 2601, Appl. Math. Comp. 217 (2011) 7238), where new class of 
^ \ soliton systems, being nonlinear integrable couplings, was introduced. The general form of 

solutions of the considered class of coupled systems is described. Moreover, the decoupling 
procedure is derived, which is also applicable to several other coupling systems from the 
literature. 



^ ■ 1 Introduction 

00 

\o 

psj ■ For a given nonlinear integrable dynamical system there usually exists many different integrable 

■ extensions which equations of motion take a triangular form. They are usually called triangular 
, systems. For Liouville integrable nonlinear ODE's the interesting class of triangular systems is 
I given in fl]. For soliton systems we know more examples. Triangular extensions of the KdV 

system was considered in [2]. In [3] there are constructed linear extensions of soliton systems, 
also taking the triangular form, being by the author called 'dark equations'. Another class of 
^ I linear extensions the so-called linear integrable couplings of soliton systems was introduced in [4j 

■ and then developed in |5| and many other papers. Recently, the theory of nonlinear integrable 
couplings of ordinary soliton systems was presented in [6] and [7| . 

In the present paper we introduce very natural triangular extensions of basic integrable 
nonlinear systems, which from the construction are also integrable. Moreover their solutions are 
uniquely determined by the solutions of the basic equation. This extensions are made on the level 
of evolution equations by a modification of the algebra of dynamical fields. The construction is 
motivated by nonlinear integrable couplings recently introduced in [6l[7]. We also propose the 
decoupling procedure for the considered class of integrable couplings including those from P, T] 

and tainiiin]. 

In Section [2] we introduce the algebra of coupled scalars, which is the underlying algebra for 
the nonlinear integrable extensions defined in Section [3l We derive soliton integrable couplings, 
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both of field and lattice type. In Section HI the general form of solution of the coupled systems is 
obtained. An example of soliton solutions of the nonlinearly coupled KdV system is presented. In 
Section [5] we derive a matrix representation of the algebra of coupled scalars and in consequence 
the matrix Lax representations for nonlinear couplings constructed in previous sections. Finally, 
in Section [6l we prove that any member of the constructed family of coupled systems separates 
into copies of the original soliton system. We also show the source of the decoupling procedure. 



2 The algebra of coupled scalars 

Consider n-dimensional vector space over field of real numbers ^ 
by the following multiplication 



We define an algebra structure 

(1) 



were are the basis vectors. Let a = Y17=i then 

f ai\ fbi\ / ci 



ab 



where 



We find that the value of the coefficient q is given by Uibi plus terms depending on lower order 
elements, ak,bk with k < i. Therefore, we call this algebra as an algebra of coupled scalars. This 
algebra is unital, commutative and associative, and are idempotent elements, what follows 
immediately from the definition ([1]). The unity element is ei. 
For n = 4 we have 



/ai\ 












b2 




as 














v 



aibi 

02^2 + 0261 + ai62 
03^3 + as(f'i + ^2) + (oi + 02)^3 



In the algebra of coupled scalars, for product of m elements the following formula holds: 

/a\\ (a^\ 



b. 



\bn, 



where 



m / k \ m /fc— 1 
j=l \r=l / j=l \r=l 



(2a) 



(2b) 



3 Nonlinear couplings of soliton systems 

Consider a commutative and associative algebra, with respect to the ordinary dot multiplication, 
of smooth functions on with m derivations -M^ : C°°(]R'") C°°(]R'"). Let us construct its 
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coupled counterpart C^(M™), that is an algebra of functions 



f(x) = /i(x)ei + ... + fn{x)er, 



where x = (xi, . . . ,Xm), taking values in the algebra of coupled scalars. So, it is commutative 
and associative algebra with respect to the multiplication ([T]) and the derivations in C^(M'") 
can be defines by 

d d 



— := i = 1, 

OXi OXi 



These derivations are well-defined as 



d 



dxi 



d 



and 



dx 



fr := — ei • f 



Here we will concentrate on the coupled extension of nonlinear integrable PDE's in (1 + 1)- 



dimension in the evolutionary form. Consider an one-field soliton system 

Ut = K[u] = K{u,Ux,Uxx,...)- 

Its extension to the system of coupled PDE's takes the form 

\it = K[u\ = K{u, u^,Uxx, ■■■), 

where 



(3) 



(4) 



u = uiei + ... +Unen 



( 



U\ = u, 



and in (jH the ordinary dot multiplication is replaced by ([T]). 
Lemma 1 The system of coupled equations ^ takes the form 

{ui)t = K [ui] 



{uk)t = K 



where K[-] is the same as in ([3]). 



.4 = 1 



K 



k-l 

,«=i 



k = 2, n, 



(5) 



The proof follows from the power series expansion of K[u] and the relation ([2|). In fact, it 
is sufficient to consider homogenous terms as 



m / k 



=1 \r=l 



m / k—1 

HE 

s=l \r=l 



m > 0, 
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where means k-ih coefficient of Uj^^. • . . . • Uj^j, and ig are arbitrary nonnegative integers. 

Note that the equation ([5]) reconstructs for ui the basic equation. Moreover, for hnear 
PDE's the procedure is trivial as all related coupled systems are just copies of the basic equation. 
The extension to the multi- field systems is obvious. 

As an instructive example let us consider the n-coupled KdV (nc-KdV): 

ut = U^xx + 6u- 



{U2)t = {U2)xxx + QU2{u2)x + HuiU2)x 



(6) 



n-1 



{Un)t = {Un)xxx + ^Un{Un)x + ^^^iUkUr. 



k=l 

It is a particular case of the triangular systems. Obviously in dispersionless limit we get n- 
coupled dispersionless KdV (nc-dKdV). 

The presented construction is also valid for differential-difference systems. If, for instance, 
instead of the derivative one defines on C^(M'") a shift operator T^e- — I ^^r^ 9 

then 

/ fi{xi,...,Xk + l,...,Xm) 

Tkei • f = : 

As an example we consider nc-Voltera differential-difference equation: 
v(x)t = v^x) ■ [v(x + 1) — v(x — 1)] 
t 

vi{x)t = vi{x)[vi{x + 1) - vi{x - 1)] 

V2{x)t = V2{x)[v2{x + 1) - V2{x - 1)] + V2{x)[vi{x + 1) - Vl{x - 1)] 
+ Vl{x)[v2{x + I) - V2{x - 1)] 

(7) 



Vn{x)t = Vnix)[Vn{x + 1) - Vn{x - 1)] + Vn{x) 
+ (^Yl '^»(^)^ + 1) - Vn{x - 1)]. 



'n-1 



n-1 



Vi{x + 1) - ^ Vi{x - 1) 



1=1 



The two field case, that is 2c-Voltera, is equivalent to the lattice integrable coupling system 
from [6]. 

Each coupled soliton system have a Lax and zero-curvature representation, that can be 
obtained by replacing an underlaying algebra by its tensor product with the algebra of coupled 
scalars. For the above examples, Lj = [L, A], where the right-hand side is the usual commutator, 
are equations in the algebras of pseudo-differential and shift operators, respectively, over algebra 
of coupled functions. Actually, for the nc-KdV ([6]) the Lax pair is given by 



d^ei + u, A = 49^ei + 6u • dxGi + Su^; 



(8) 
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and for nc-Voltera ([7]) by 

L = Tei + v(x)-T-^ei, A = T^ei + v(x + 1) + v(3;). (9) 

As an example of coupled extension of two-field system let us consider 2-coupled AKNS 
(2C-AKNS), n = 2: 

1 1 

Pt = -^Pxx + p • p • q, qt = ^<ixx - p • q • q 

t 

1 2 1 2 

Pt = -^Pxx +P q, qt = ^qxx - pq (10) 

1 o 9 9 

rt = —-Txx + p s + 2pqr + 2prs + qr + r s 
St = ^Sxx - ps^ - 2pqs - 2qrs - q^r - rs^, 



where 



:p\ 



In fact the system (|10p was obtained in [8j in a different way and is equivalent to the one from [7] . 
Its zero-curvature representation will be presented in Section [H 

4 Solutions of coupling systems 

The solutions of the coupled systems (jl]) are completely determined by solutions of the basic 
equations. 

Theorem 2 Assume that S^,...,S^ are arbitrary different solutions of the basic equation ut = 
K[u\. Then, the solution of coupled system ^ is given by 

ui = S^, Uk = S^-S''~\ k = 2,...,n. (11) 

Moreover, any solution of the coupled system is of the form ill]) . 

The first part of the theorem is evident, after plugging (lllh to ([5]) we get 

The proof of the second part is by induction. If the solution of the basic equation is ui = S^, 
then the first coupled equation is 

{U2)t = K[U2 + S^] - K[S^] ^ {U2 + S^)t = K[U2 + S% 

thus U2 + = S"^, where S'^ is another solution of the basic equation. The rest of the induction 
is obvious. The consequence of the theorem is that for coupled PDE's any n different solutions 
of the basic equation build up an appropriate solution of the coupled system. 

Let us illustrate the whole construction on the example of solitons of nc-KdV ([6]) in the 
Hirota form. Let us start from 1-soliton solution. Consider n one-soliton solutions of the KdV 
with the same wave velocities and different phases 7^: 

= 2dx log(l -|- /■'), f-' = exp(kx + k^t + 7j), k, jj = const, j = 1, n. 
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Figure 1: Coupled 1-soliton, Ui = — ^ / !)• 
Then, 1-soliton solution of (j6]) is given in the form 

m = S^ = 2dl log(l + /I), = - S^-' = 2dl log [ l^jl^ > i = 2, n. 

Obviously, u\ is the ordinary KdV soliton, while remaining Ui are coupled solitons (illustrated 
in Fig. [1]), that differ among themselves by phase arguments. 
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Figure 2: Coupled 2-soliton, m = - S'-^ (i / 1). 
The coupled 2-soliton solution of nc-KdV ([6]) is of the form 
ui = S' = 2dl log(l + l\ + fl + aflfl), 

i+f(+n+af(n \ 



Uj = - S^-^ = 2dl loi 



where 



fl = exp(/cix + kft + 7ij), = ex.p{k2X + kit + 72^) 



3 = 2,...,n, 
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and a = {ki — ^2)^/(^1 + ^2)^- Notice that ui is the ordinary 2-sohton solution of the KdV, 
while remaining Ui are coupled 2-soliton solutions, which interaction is presented in Fig. [2j 



5 Matrix representation of the algebra of coupled scalars 

Consider n quadratic matrices E^. of dimension n x n: 



1 if j = k, i < k 
1 if j = i, i > k 
otherwise, 



where k = 1, . . . ,n. The matrices Ej constitute generating elements of commutative and asso- 
ciative sub-algebra of triangular matrices. Let us call the matrix algebra spanned by E^ as a 
pseudo-scalar algebra of matrices, ps{n). 
For example if n = 4: 



El 



E3 

Then, a typical element for n = 4 is 

n 

A = ^ OjE, 



/I 












/o 


1 










1 








E2 — 





1 














1 











1 





\o 








V 












1/ 







1 














1\ 








1 





E4 = 











1 








1 














1 


Vo 








V 












1/ 



/ai 


02 


as 


04 \ 





ai + 02 


03 


CI4 








ai + 02 + as 


04 


Vo 








ai + 02 + as + 04/ 



(12) 



Lemma 3 The algebra ps{n) is a matrix representation of the algebra of coupled scalars defined 
by the multiplication ([T]). 

The lemma follows immediately by showing that 

Tp Tp 

^i^j — ^j^i — ■CJmax(i,j) 

for i, j = 1, . . . ,n. 

The Lax representation L^ = [L, A] of nc-KdV ([8]) and nc-Voltera ([9]) in the matrix algebra 
ps{n) are respectively given by 



where 
and 



L = a^Ei + C/, A = 4a^Ei + 6C/a^.Ei + 3f/^, 

U = uiEi + n2E2 + ... + tinE„, 
L = TEi + V{x)T-^Ei, A = T^Ei + V{x + 1) + V{x), 
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where 

V{x) = Vi{x)Bi + V2{x)E2 + ... + Vn{x)Bn. 

Taking the tensor product of some Lie algebras, like in the above examples of pseudo- 
differential and shift operators, with the pseudo-scalar algebra ps{n) one can derive several 
examples of the coupled extensions of known integrable systems. 

This is in fact the case of the coupled AKNS systems (fTUj) . where its zero-curvature equation, 

Lt- W, + [L,W] =0, (13) 

is from the Lie algebra being tensor product of ps(2) with the loop algebra of sl{2). 
The respective generating operators are given by 



where 



and 



where 



xr ^ Vs 



W = El O PTi + E2 O W2 



Wi W2 
W1 + W2 



e 



with 



> qX+^Qx X^-^pqJ \f -d 

d = a X^ + -{1 — a){ps + qr + rs) — -apq, 
e = ((1 - a)r - ap)X + ^ap^^. - ^(1 - a) r^, 

/ = ((1 - a)s - aq)X - ^^Qx + ^(1 - a) s^. 

We have left here some freedom related to the parameter a being integration constant from the 
computation of W. For various details on derivations of integrable systems we send the reader 
to [TT] and references therein. Then from the zero-curvature equation (|13|) we get the following 
system 

1 2 1 2 

Pt = --Pxx +p q, qt = -qxx - pq , 

^ I /I ^ \ I ^ ^ I o ^ I o^„„ I 2 I ^2„ I „2, 



rt = -apxx - ap q+{l- a) y-- r2x + 2pqr + 2prs + qr +p s + r sj (15) 

St = -^afc + apq'^ + {a - 1) ^-^ S2x + 2pqs + 2qrs + q^r + ps^ + rs^ 

The 2c- AKNS (jiup (also obtained in [S]) one obtains for a = 0. Whereas, the coupled AKNS 
obtained in [7j is generated for q = — 1. This system is in fact a linear composition of the 
symmetry (jlOh (the case of a = 0) and the one obtained from (llSp for a = 1. 

In the papers [6]-[10j there are also other examples of integrable couplings for field and 
lattice soliton systems of the same class. All of them are generated by the spectral Lax operators, 
from the Lie algebra ps{2) (8) s/(2)[A, A~^], being in the form ()14p . however with different entries 
of Ui and 1/2- 
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6 Decoupling procedure 

Consider the transformation of the pseudo-scalar algebra ps(n) in the form of a similarity relation 

r(A) := S-^AS, 

where 



and 



1 for i < j 
for i > j 



1 for j = i 
— 1 for j = i + 1 
otherwise. 



Then, this transformation is apparently an isomorphism of matrix algebras. One finds that 



1 for j = i, j > k 
otherwise. 



In particular, for n = 4 



S 



4 








/I 
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1 


1\ 
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^0 








1/ 



s 



-1 



/I 


-1 










1 


-1 











1 


-1 










1/ 



and for A given by p2] 
T(A) 



/ai 

ai + 02 

VO 



\ 



01 + 02 + 03 

01 + 02 + 03 + 04/ 



The above transformation naturally extends to any tensor product of some algebra q with 
ps{n), that is ps{n) (g) g. Let 



(16) 



k=l 



where belong to g. Then, the transformation extends to elements of the form ()16p by the 
formula 



r(L) :=^r(Efe)0Afc. 



(17) 



k=l 



Since this transformation is an algebra isomorphism it preserves the Lie algebraic construction 
of the coupled systems. Combining ()17p with the transformation of components: 



k = l. 



(18) 



1=1 
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we decouple the whole construction into n independent copies, that is ps{n) (8) decouples into 
the direct product of n copies of q. 

As consequence of the above decoupling procedure, the linear transformation to new field 
variables 

k 

Uk = ^Ui k = l,...,n, 

i=l 

separates the coupled equations ([l])-© to n copies of the basic equation ([3]) in new fields Uk, 
that is 

{uk)t = K[uk] k = l,...,n. 

In particular this is the case of nc-KdV ([6]) and for nc-Voltera ([7]) systems. 

In the case of the coupled AKNS system (jlOp or (jlSp generated by the Lax operator (jl4p 
the decoupling transformation (jlSp takes the form 





= Ui 


Pi 


= P 






=^ P2 


= p + r 




= Ui + U2 




= q + s. 




92 



Qi = q 



As result the coupled AKNS separates into two copies of the standard AKNS system. 

All coupled triangular systems considered in previous sections as well as these from [6]- 
[10] (as they construction is based on the Lie algebra ps{2) s/(2)|A, A~^] separate (decouple) 
into copies of the basic equations in the same manner. In fact the same situation is in the 
general case, when the underlying Lie algebra has the form ps{n) Q. This follows from the 
fact that the above transformation is in the form of similarity relation which naturally preserves 
Lie commutator, trace form. Lax and zero-curvature representations, Hamiltonian structure and 
other related structures. 
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